Slow plasmon modes in polymeric salt solutions 
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The dynamics of polymeric salt solutions are presented. The salt consists of chains A and B, 
which are chemically different and interact with a Flory-interaction parameter x> the A chain ends 
carry a positive charge whereas the B chain ends are modified by negative charges. The static struc- 
ture factor shows a peak corresponding to a micro phase separation. At low momentum transfer, 
the interdiffusion mode is driven by electrostatics and is of the plasmon-type, but with an unusu- 
ally low frequency, easily accessible by experiments. This is due to the polymer connectivity that 
introduces high friction and amplifies the charge scattering thus allowing for low charge densities. 
The interdiffusion mode shows a minimum (critical slowing down) at finite k when the interaction 
parameter increases we find then a low k frequency quasi-plateau. 



Scattering methods are a powerful tool to study the structure and phase behavior of multicomponent polymer 
mixtures and were the subject of intensive investigations [Q. In general, polymer blends are characterized by a low 
entropy of mixing and relatively high enthalpic interactions and polymers are thus to a large extend immiscible 
iBoJliul- ^ n effective way to enhance the compatibility of polymer blends, consists to introduce long range Coulomb 
forces between particles by adding electrostatic charges . 

Quasi elastic scattering experiments allow, in principle, to measure the decay of various concentration fluctuations 
in the solution Q. The fastest relaxation mode is usually the plasmon mode related to the relaxation of charge 
fluctuations. In the macroscopic limit (k — > 0) for small ions, the plasmon frequency is proportional to the inverse 
friction coefficient £, and to the Debye-Hiickel screening parameter re 2 , i.e., T cx ksTn^C 1 ■ Although it is possible 
to detect the plasmon by spin echo techniques it remains very difficult to measure, because of low contrast. Then 
other processes such as water binding - unbinding perturbate a clear detection of the plasmon. To find experimentally 
accessible plasmon modes in reasonable time scales larger friction coefficients, low charge densities, and higher contrast 
are required. We suggest to use polymeric salts, where anions and cations are linked to the end of two chemically 
different polymeric backbones. To be more precise, chains A and B interact by a Flory - \ parameter, which allows for 
composition fluctuations Many hydrophilic polymers involve hydrogen bonding and may exhibit a more complex 
behavior. These effects can be limited by an appropriate choice of temperature and polymers, e.g., end- modified 
POE at room temperature. The chain friction is a power of N ||] larger than the monomer friction of the order of £ 
introduced above. 




FIG. 1. Polycation A and polyanion B in solution. 



To benefit from the polymeric effects small ions should be absent, which is possible in a one to one A-B mixture, 
this point will be discussed later. Electroneutrality forbids macrophase separation as x increases and microphase 
separation at a finite wave vector k* is anticipated. Large scale composition fluctuations violate electroneutrality and 
their relaxation is primarily driven by electrostatics. We will indeed show that at low wave vector the relaxation of 
composition fluctuations is of the plasmon type but at an unusually slow time scale. Though critical slowing down 
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closer to the microphase separation is interesting, we will mainly concentrate on systems far from the transition 
where excluded volume effects do not interfere with the plasmon effects. For the weakly charged chains we have 
in mind that the value for k is small. On the other hand the polymer concentration is high and high contrast in 
scattering experiments can be achieved. Symmetric polymer salt solutions with moderately incompatible backbones 
(e.g. hydrogenated and deuterated) are thus good candidates to detect slow plasmon modes. 

The static properties of the polymeric salts in equilibrium are well described by the random phase approximation, 
since we assume a sufficiently large overall concentration. The chains are assumed to be of Gaussian connectivity, the 
monomers interact with an excluded volume interaction of a strength given by Vij , and the electrostatic interaction 
between the ionic groups. In the latter the strength of the Coulomb potential is expressed in term of the Bjerrum 
length: Ib = f3e 2 /Aire, where e and e are respectively the electric unit charge and the dielectric constant. The different 
chemical nature of the chains is commonly expressed by a Flory \ - parameter, Vab — V + X- F° r the purpose here it 
turns out convenient to use linear combinations of the collective variables to study the usual plasmon (total charges) 
modes and the interdiffusion (relative A and B fluctuations). It is straightforward to calculate the corresponding static 
structure factor matrix S(k) which is expressed as a function of the bare structure matrix So(k), corresponding to 
the single chain conformations in the mixture, and the interaction matrix U; n t(k) in the form 

S(k)" 1 =S (k)- 1 +U int (k). (I) 

The interaction matrix U; n t(k) contains short (excluded volume) and long (electrostatic) range interactions. The 
elements of the first matrix are the usual parameters Vij . The relevant structure factors are S zz (k) (charge / charge 
structure factor), S' II (k) (ionic strength / ionic strength structure factor), S^(k) (concentration / concentration 
structure factor), S* xx (k) (composition / composition structure factor), S l ^(k) (ionic strength / concentration structure 
factor), 5 zx (k) (charge / composition structure factor). In the case of the symmetric salt where chains A are charged 
positively and B charged negatively, 5* xx (k), S zz (k) and 5 zx (k) have a peak which increases as the value of the Flory 
parameter x increases. In this case we find a microphase instability with a characteristic finite wavevector k*. The 
microphase separation at non- vanishing vector k* occurs for a sufficiently strong incompatibility. For x — or \ finite 
but sufficiently small, entropy effects are dominant and favor mixing. 

In light or neutron scattering experiments, one can measure directly the matrix of dynamical functions Sij(k, t). 
To compute the relevant elements of the dynamic structure factor matrix we start from the generalized Langevin 
equation, which describes the (short) time evolution of S(k, t) and, neglecting memory effects, lead us to a single 
exponential decay of S(k,t), S(k, t) — e _n ( k )'S(k). The decay rate, which is experimentally available is given by the 
first cumulant matrix f2(k) which is related to different correlation functions. Hydrodynamic interactions between 
monomers are taken into account at the Oseen level [[) f . 



«<k) ,V,3 /^S(k + q) fc2 M s-'(k) (2) 



The eigenvalues of fl(k) lead us to the decay rates of thermal fluctuations in the solution. In a first step the eigenvalues 
of the frequency matrix eq. @) have been computed numerically. 

Let us first concentrate on the fast (usual) plasmon mode, which corresponds to the relaxation of charge fluctuations. 
At high wavevector limit, the relaxation is diffusive as for a monomeric salt and the electrostatic interactions are 
irrelevant. In the low k limit, the symmetric salt behaves as an usual monomeric salt. The relaxation of charge 
fluctuations of size larger than k _1 is driven by electrostatics. Since the chain ends are different from the rest of the 
chain monomers, an additional relaxation mode is expected. At high wave vectors, the relaxation is diffusive and at 
larger length scales, the relaxation is rather driven by chain elasticity. The slowest mode corresponds to the relaxation 
of monomer concentration fluctuations and is similar to the mode observed in binary polymer-solvent solutions. This 
mode is independent of the strength of the electrostatic interaction and, depends mainly on the excluded volume V. 
It is insensitive to the Flory parameter x- At small length scales, the relevant friction is proportional to the probed 
length scale 1/fc and in the limit of the large length scales, the relevant friction is proportional to the correlation 
length £. 

Let us now turn to the interdiffusion mode. If the charges were absent the interdiffusion would show the classical 
critical slowing down of the composition fluctuations close to the macro phase separation. Charges stabilize the solution 
against macrophase separation and the relaxation of composition fluctuations at low k is dominated by electrostatics: 
the relaxation frequency has a finite value in the macroscopic limit. This behavior is characteristic for plasmon modes. 
As can be seen on the Fig. 2 below the frequency is very low (of the order of I0 4 Hz, for the typical values shown). 
At wavevector zero a weak (not visible at the scales of Fig. 2) splitting occurs, due to the dependence of the mobility 
on the Flory- x parameter. The behavior found here is somewhat similar to block copolymers. The usually observed 
critical slowing down is shifted to the local scale k* . Thus the corresponding relaxation frequency shows a minimum 
for high enough values (higher than chosen in the figure) of the interaction parameter x- 
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FIG. 2. Variation of Intcrdiffusive mode versus the wave vector k. Dashed curves correspond to 
system (2X2) Continuous curves correspond to system (4x4) for different values of \ ( 0, 0.1 and 
0.24 ), N = 10 4 ,6 = l B = 710" 8 cm, <f>a 3 /N = 3.6710" 6 . 

The analytical calculation of the full 4x4 system involves tedious algebra and produces long and non transparent 
formulae. In the low wavevector limit the precise position of the charges is not of specific interest PL Therefore it is 
convenient to smear out the charges. Thus it is possible to replace the system described above by an effective two chain 
model with a charge / = 2/N per monomer. The static composition structure factor S^x can then be represented as 

s =VL{ i \ (3) 

2 \(2-x0AO + fc 2 i?2 + 4K 2 A 2 / ' 1 ' 

where k 2 = 4ttIb(2c/))/N. If the chains are uncharged, the classical case of a binary polymer blend is recovered. 
The term proportional to 1/fc 2 stems from the bare Coulomb interaction between the charges. 5 XX diverges then at 
k = and x<t>N = 2. The charges themselves enhance compatibility and yield a microphase separation transition at 
X(j)N = 2 + AnR g . The structure factor S xx shows a characteristic peak at a finite wavevector k* where k* 2 — 2n/R g . 
For understanding the dynamical results obtained by numerical computation of the full four component system it is 
further convenient to rewrite the structure factor to 

cj)N i r i i 
~ 1 1 + k 2 d ~ ~i + k 2 s; 2 



S KX - — -j= \ - ; u2t2 - - ■ 2 \ (4) 



where we have defined the discriminant A = (2 — %0iV) 2 — (4K_R g ) 2 . In this representation two important length scales 
show up, i.e., £ 2 2 = (2i? 2 )/(2 — x4>N ± V~A). It is immediately seen that £i£ 2 = k*~ 2 . For weak incompatibility, i.e., 
X<pN < 2 — 4/ci?g, the correlation function in real space consists of two exponential decaying functions exp(— r/^i^), 
see eq (|]). Closer to the microphase instability, 2 — 4nR g < x4>N < 2 + 4K_R g , the real space structure factor contains 
decaying oscillatory parts. In the latter regime the plasmon behavior of the composition mode at low k will turn out 
to be hidden by the composition fluctuations. Moreover renormalization theories beyond the RPA are required. 

The dynamics can be understood within the same framework. The relaxation frequencies can be computed along 
the lines discussed above. It is sufficient to understand the behavior of the relaxation frequencies from figure (1) at 
low values for the wavevector. Most interesting is the limit at zero wavevectors 

fi xx (k = 0) = (I + I)" 1 , (5) 

which has to be discussed in several limits. The first is the limit of x = 0, i.e., no composition fluctuations and 
therefore no coupling between thermodynamic and electrostatic interactions. For this case we find the simple result 

k B T 4k 2 



1 This statement is less obvious for the mobility integral occurring in eq.(^|). However it can be checked afterwards that the 
integral is not dominated by the large k range. 
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where kR 9 << 1 was taken into account. The physical interpretation here is indeed simple: The relaxation frequency 
corresponds to a plasmon mode of <fi/N macroions carrying two elementary charges each, whose hydrodynamic inter- 
action is of the range of the radius of gyration R g (corresponding to the Oseen - tensor 1/r). Composition fluctuations 
do not couple in this limit, thus the non appearance of any wavevector k* is natural, in contrast to the upper boundary 
X^iV — 2 — 4nR g . There the first sign of the micro structure shows up in terms of k* . Indeed we find 

k T 

CWx0^ = 2-4 K i? g ) = -^-(fc*) 3 /2 (7) 

In this regime of the interaction parameter the dynamics is ruled by the intrinsic length scale k* of the microphase 
separation transition. The plasmon character of the composition mode is reflected by eq. @. For its observation a 
weak incompatibility like between deuterated and hydrogenated polymers is necessary. In the other limit which is 
given by eq.(0) the plasmon character is already blurred by the microphase separation. The slow plasmon mode is 
predicted in a range accessible by quasi elastic light scattering and spin echo techniques. The dynamical study has 
been limited to short time scales where the decays are exponential. Entanglements have been ignored completely. With 
respect to our studies, they will only further slow down all polymeric modes. The measurements are carried out in 
aqueous solution, and small amount of salt and ions could destroy this behavior. We carried out a model calculation, 
which involves these ions as individual component in the Langevin dynamics. The small ions couple to the electric 
charges bound at the chains, via a Coulomb potential. The dynamics of these ions with respect to all polymer modes 
involving the composition and concentration is a power of N faster, N being the chain length. Therefore the fast small 
ions produce a screened electrostatic potential of the Debye - Hiickel type, where the screening constant is determined 
by the concentration of the remaining small ions. The relevant concentration in water is often of order 1CP 6 , slightly 
higher than imposed by self-dissociation, and the corresponding wave vector is about k = I0 4 cm -1 . Thus we expect 
that at larger wave vectors the slow plasmon mode can be observed. 

Acknowledgement: The funding of this work by the Deutsche Forschungsgemeinschaft, Max-Planck-Society and 
L.E.A. are most gratefully acknowledged. We thank Henri Benoit, Michel Rawiso, Jean-Frangois Joanny, Mustapha 
Benmouna and Burkhard Diinweg for many helpful discussions. 



[1] J. S. Higgins, H. Benoit, Polymers and neutron scattering, (Oxford University Press, Oxford, 1994). 
[2] T. Ohta and K. Kawasaki, Macromol. 19, 2621 (1986). 
[3] F. S. Bates, Science 251, 898 (1991). 

[4] T.A. Vilgis and J. Noolandi, Macromolecules 23, 2941, (1990) 
[5] G. H. Fredrickson, Macromol. 26, 2825 (1993). 
[6] M. G. Brereton, T. A. Vilgis, Macromolecules 23, 2044, (1990). 
[7] T. A. Vilgis and P. Haronska, Macromol. 27, 6465 (1994). 

[8] P.G. de Gennes, Scaling concepts in polymer physics, (Cornell Unviersity Press, Ithaca, 1979). 
[9] M. Doi and S. F. Edwards, The theory of polymer dynamics, (Oxford University Press, Oxford, 1986). 
[10] T. A. Vilgis, M. Benmouna, J. Chem. Phys. 95, 1276 (1991) 



4 



